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Abstract: A transformation is obtained which completes the unification of 

quadrirational Yang-Baxter maps and known intcgrablc multi-quadratic quad 
equations. A generalised discrete dynamics emerges by combining theory from 
these two classes of quad-graph models, it is encoded in a birational representa- 
tion of a particular sequence of Coxeter groups. In this setting the usual quad- 
graph is associated with a subgroup of type i?C„, and is part of a larger and 
more symmetric ambient space. The model also defines, for instance, dynam- 
ics on a triangle-graph associated with a subgroup of type An, as wejl as finite 
degree-of- freedom dynamics, in the simplest cases associated with and Eg 
affine subgroups. Underlying this structure is a class of biquadratic polynomials, 
that we call idempotent, which express the trisection of elliptic function periods 
algebraically via the addition law. 

1. Introduction 

The point of departure in this paper is the investigation of integrable quad 
equations. It will be demonstrated that the primary model in the multi-quadratic 
class that was discovered recently by the author in collaboration with Nieszporski 
[1] , is related in a natural way to the primary model of the quadrirational Yang- 
Baxter maps that have been introduced earlier by Adler, Bobenko and Suris [2]. 
The existence of this relationship clarifies the position of the new model with 
respect to known quad-graph models and transformation theory in [3-15]. The 
subsequent focus of this article is on a particular repercussion of this connection, 
namely a new interpretation in terms of the idempotent class of biquadratic 
polynomials. 

The idempotent biquadratic has its roots in the theory of elliptic functions; it 

is the formula for trisection of the periods expressed algebraically via the addition 
law. Interesting properties of this formula in part correspond to integrability of 
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the aforementioned quad-graph models, but go also beyond that, in particular 
leading to a natural generalisation of the quad-graph dynamics. 

A practical way to understand the additional structure is in terms of the braid 
group representation emerging from the Yang-Baxter maps [16]. In the quadri- 
rational case the braid-type generators are self-inverse, and the representation 
is only of the trivial part of the pure braid group, i.e., the symmetric group. 
But here we identify this as a subgroup of a larger ambient group. The precise 
characterisation of this larger group is one of the main technical achievements 
of this paper, it is a Coxeter group with unexpectedly rich structure (cf. Figure 
10). For instance, in the sequence of groups formed by an increasing number of 
generators, the last finite case is the exceptional finite reflection group Eq, and 
the associated lattice geometry coincides with the configuration of the 27 lines on 
a cubic surface. The general group structure bears closest resemblance to groups 
whose birational representation forms part of the basic framework in which the 
Painleve equations and their generalisations may be understood [17 -22]. 

We proceed as follows. Section 2 introduces the aforementioned quad-graph 
models and establishes the new relationship between them. The idempotent 
class of biquadratics is introduced in Section 3 with additional material given in 
Appendix A. Section 4 connects the quad-graph models with the idempotents via 
the notion of quadrirationality. A natural closure property of the idempotents, 
and an associated integrable triangle-graph dynamics, is established in Section 
5. The nature of the relationship between the quad-graph and triangle-graph 
dynamics is a quite delicate matter, motivation behind a method to unify them 
is given by a local consideration of the problem in Section 6, and a constructive 
approach which completes the solution globally, based on vertex maps, is given 
in Section 7. The solution emerges in the form of the aforementioned birational 
group structure, which is interpreted from the point of view of lattice-geometry 
and dynamics in Section 8. Some consequences of the results here, and natural 
questions which emerge, are discussed in Section 9. 



2. Preliminciries, Fj and Q4* 

The two relevant classes of models are the quadrirational Yang-Baxter maps 
and integrable multi-quadratic quad-equations. We recall the primary model 
from each class and connect them by a new Miura-type transformation. 

Rational mappings with the Yang-Baxter property were constructed by Adler, 
Bobenko and Suris [2] on the basis of a simpler property termed quadrirational- 
ity. The primary model obtained was denoted Fj: 

(1 - (3)u + P - a- (1 - a)v 
(1 - a)/3w + {a- p)uv - (1 - I3)av ' 
, ^ {l-a)v + a- p)u ^ ' ' 

y — [in 

^ {1- P)av + {P-a)vu-{l-a)Pu' 

It is natural to view such models as dynamical systems on a quad-graph [23— 
25], see Figure 1(a); variables assigned to edges are governed by the system 
(2.1) on each quad, whilst the essential parameters of the model, a and 13, are 

associated with characteristics. This is unambiguous due to invariance of (2.1) 
under permutations u ^ u, v ^ v and {u,u,j3) O {v,v,a), or in other words 
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Fig. 1. (a) Planar quad-graph, (b) quad with variables on edges, (c) quad with variables on 
vertices. Parameters are associated with the dashed lines, which represent characteristics. 

the system respects the symmetry of the quad, see Figure 1(b). The rationahty 
of the mapping plus the first two of these permutation symmetries is the essence 
of the quadrirationahty. 

A class of transformations connecting models with variables on edges, such 
as (2.1), with polynomial quad equations in which variables are on vertices, 
were introduced by Papagcorgiou, Tongas and Veselov [10] and developed by 
Kassotakis and Nieszporski [12]. It is a new transformation in this class which is 
the point of departure in this article, it involves equations on edges of the graph 
which around a single quad take the form 



where B is the following polynomial 
B{u,x,x,a) := 

(1 + cx){l + cx)u'^ - [(1 + c^)(l + xx)a + 2c{x + x)]u +{c + x){c + x)a, 



and c e C\ {0, 1, —1, i, —i} is a constant parameter. This system connects model 
(2.1) with the quad equation 



B{u,x,x,a) = 0, B{v,x,x, p) = 0, 
B{u,x,x,a) = 0, B{v,x,x, = 0, 



(2.2) 



(2.3) 




(2.4) 



where 



p=m{a), g = m(/3) 



m := a 



2c^ - (1 + c^)a 
c(l + c2)a-2c' 



(2.5) 
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Fig. 2. The Yang-Baxter property, global transfer of data along characteristics is independent 
of the order in which they cross. This is equivalent to the local condition that equations for 
variables on white nodes in terms of variables on blax;k ones are the same for both of these 
quad graphs. 

see Figure 1(c). The quad equation (2.4), denoted Q4*, is quadratic in the depen- 
dent variable, but the basis ol its original construction [1] guarantees existence 
of a rational reformulation. In fact the system on edges (2.2) is a way to intro- 
duce auxiliary edge variables, the (autonomous) rational reformulation in terms 
of these variables being exactly model (2.1). Typically such reformulation would 
involve both edge and vertex variables, so relations (2.2) are especially natural. 
It means system (2.2) is also a Miura-type transformation obtaining solutions of 
(2.4) from solutions of (2.1). 

In fact these features can be explained in terms of the integrability of (2.1), 
namely the Yang-Baxter property, or cubic consistency, see Figure 2. The main 
observation is that 

B{u,m{w),m{w),a) = <^ 



thus (2.2) is really nothing but the same underlying model (2.1). The system of 

all equations (2.1), (2.2), (2.4) is associated with the non-planar cube that, due to 
the consistency property, may be formed as the union of the two graphs in Figure 
2 without breaking the well-posedness of the indicated initial value problem. Tliis 
cube is shown in projected form in Figure 3 with variables and parameters added 
explicitly. It can be verifed by calculation that as a consequence of imposing 
(2.1) on each face, the variables on the four edges through which the closed 7 
characteristic passes are related by (2.4) via the substitution 



(Note that the generic form of (2.4) will be given later in this article, Mobius 
transformation m here and in (2.5) is characterised by its action m : (0, 1, 00) 1— )• 
(— c, 1/c, — 1/c) ). Up to the same change of variables, the equations (2.2) are as- 
sociated with the four quads through which the 7 characteristic passes. Upon this 
observation the statement that (2.2) constitutes a Miura-type transformation be- 
tween (2.1) and (2.4) can be identified with the statement of the Yang-Baxter 
property of (2.1) itself. 

The theory of the simpler rational model (2.1) therefore encompasses the 
transformation (2.2) and the multi-quadratic model (2.4). In hindsight one can 




(2.6) 



X = rr){'w), X = m(tl)), x = m(w), x = m['w), 
p = m(a), q = m(/3), c = m(7). 



(2.7) 



Idempotent biquadratics, Yang-Baxter maps and birational Coxeter groups 



5 





Fig. 3. A quad-graph domain for (2.1). The four variables along the closed characteristic 7 
are related by the multi-quadratic quad equation (2.4) modulo the change of variables (2.7). 

notice that the basis of construction of (2.1) and (2.4) are in fact similar. In the 
first instance it is the preservation of rationality under permutation of the vari- 
ables of the model, and in the second it is existence of a rational reformulation. 

However, there is one particular element from the theory of (2.4) which turns 
out to give new insight into (2.1). It is known [1] that the interpretation of B as 
Backlund transformation for (2.4) implies the following discriminant property, 

A[B{u,m(w),m{w),a),u\ — 

[(1 — j)~^{w — j){w — j){a — 7) + 7(ui 4- w -I- a — 7)]^ — Ajwwa = 0, (2.8) 

where the polynomial appearing in (2.8) is the Weierstrass-typc symmetric tri- 
quadratic with discriminant polynomial r{w) — ■w{w — l){w — 7). This char- 
acterises B as being amongst a particular class of polynomials that we call 
idempotent, as shown in Appendix A. Due to (2.6), the model (2.1) is also con- 
nected with this class of polynomials, and it is this connection which is the key 
observation we develop in this paper. The idempotent class has some remarkable 
features of independent interest, and the connection with (2.1) can be established 
with greater directness after some of those features are understood. 

3. Idempotent biquadratic correspondences 

We define and give basic properties of the relevant class of algebraic correspon- 
dences. 

Symmetric biqudaratic polynomials, 

h{x, y) = co + ci{x + y) + C2xy + 03(0;^ + y^) -|- C4xy{x + y) + c^x^y"^, (3.1) 

co,...,C5 € C, provided they don't decompose into a product involving one- 
variable polynomial factors, define dynamics on orbits, that is sequences of values 

. . . ,a;_2,a;_i,a;o,a;i,a;2, . . . e C U {00} (3.2) 

such that 

h{Xn,y)=0 <^ y e {Xn-l,Xn+l}, n S Z. (3.3) 
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The two- valued dynamics defined by the biquadratic manifests only in the choice 
of direction along the orbit. The significance of the orbits, which are sequences 
defined uniquely up to orientation, is that they may be taken as a definition of 
elliptic functions in the discrete setting. 

The self-composition of the correspondence defined by h yields another sym- 
metric biquadratic polynomial, g, via the formula 

Tes2[h{x,y),h{y,z),y] = {x- zfg{x,z), (3.4) 

where res2 denotes the resultant of quadratic polynomials^. The other factor 

appearing in (3.4), {x — z)^, is due to the fact that h is symmetric. This notion 
of self-composition allows to define the idcmpotcnt class. 

Definition 1. We refer to h (3.1) as idempotent if its self-composition g, defined 
by (3.4), is a non-zero scalar multiple of h. 

Note that if h has a one-variable polynomial factor, or it drops in degree, then its 

self-composition vanishes identically. Such degenerate cases are excluded in Def- 
inition 1 by requiring the scalar multiple be non-zero. In particular, idempotent 
biquadratic polynomials have well-defined orbits. 

The utility to be found in biquadratic polynomials from this class begins with 
the following list of complementary ways to characterise them. 

Proposition 1. Assume the symmetric biquadratic polynomial h in (3.1) is not 
a scalar multiple of {x — y)"^, then the following are equivalent, 

(i) h is idempotent, 

(a)* all orbits of h are three-periodic (x^+z = x„ in (3.2)), 
(Hi)* there exists a non-constant three-periodic orbit of h, 

(iv) * the coefficients of h satisfy the constraint c\Ci + C3 = C0C5 + C2C3, 

(v) * h satisfies the four-variable polynomial identity 

{w - z)h{w, z){x - y)h{x, y) + {w - y)h{w, y){z - x)h{z, x) 

-\- {w - x)h{w,x){y - z)h{y,z) =0, (3.5) 

(vi) h can be written in the form 

h{x,y) = ''^^^''^y^-''^''^''^y\ (3.6) 
x-y 

where ri , r2 are cubic polynomials with no common roots. 

Here the * indicates addition of the non- degeneracy condition on h that its self- 
composition, g in (3.4), does not vanish identically. 

Proof. Deleting every other point from an orbit of h yields an orbit of its self- 
composition g defined by (3.4), considering this and the fact that orbits of g and 
h coincide due to (i), leads to (ii)*. 

Existence of a non-constant orbit is a consequence of the assumption that h 
is not a scalar multiple of (x — y)^, therefore (ii)* implies (iii)*. 

^ More precisely res2 Ip{x), q(x), x] is the determinant of the Sylvester matrix of p and q with 
the expUeit assumption that both p and q arc degree two polynomials, allowing a possibly zero 
coefficient of the second-degree monomial term. 
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Fig. 4. An orbit of an idempotent biquadratic: h{xi,y) = <^ y E {xj,Xh}, {hj,k} = 
{0,1,2}. 

To see that (iii)* implies (iv)* denote the elements of the three-periodic orbit 

mentioned in (iii)* by xq, xi and X2, at least two of which are distinct, and 
consider the resulting three equalities between polynomials in y: 

Hxo,y) = {c3 + C4X0 + C5Xo)(y - xi){y - X2), 

Hxi.y) = (c3 + C4X1 + czxl){y - X2){y - xq), (3.7) 

Hx2,y) = (c3 + C4X2 + c^xl){y - xo)iy - xi). 

Substituting the explicit form of h (3.1) followed by elimination of xq, xi, X2 
yields the condition on the cocfRcicnts which appears in (iv)*. 

Substitution of (3.1) directly shows that polynomial identity (3.5) is satisfied 
as a consequence of the condition on the coefficients appearing in (iv)*, so (iv)* 
implies (v)*. 

Choosing w and z to be fixed constants such that 77^ := {w — z)h{w, z) 7^ it 
is clear that 

-(w — x)h{w,x) = ri{x), -{z — x)h{z,x) = r2(x), (3.8) 
V V 

are cubic polynomials in x, and writing (3.5) in terms of these one obtains (3.6). 
Furthermore ri , r2 have no common roots due to the non-degeneracy assumption 
on h, therefore (v)* implies (vi). 

Assuming h is of the form (3.6) specified in (vi), a direct calculation shows 
that 

res2[h{x,y),h{y,z),y] = {x - zf ^ih{x , z) , /i := res3[ri(y), r2(y), y]. (3.9) 

The scalar multiple ji appearing in (3.9) is non-zero if and only if ri,r2 in (3.6) 
have no common roots (resa excluding also the possibility of 00 being a common 
root, which corresponds to a drop in degree of both polynomials), therefore (vi) 

implies (i). 

Equivalences (ii)* and (iii)* connect the idempotent class of biquadratics with 
the trisection of the elliptic function periods. The three-periodic orbits will play 
the central role in what follows. They are specified completely by the un-ordered 
triplet of participating points, orbits being determined only up to orientation. 
Therefore within this class of biquadratics the un-ordered triplet itself may be 
referred to as the orbit, and, to reflect the symmetry, it is natural think of points 
in an orbit as being arranged in a triangle as in Figure 4. 

We remark that the idcmpotcnts occupy a similar special position in relation 
to the biquadratic correspondences as do Mobius involutions to the group of 
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Mobius transformations. This was observed in [26] based on the similarity of 
characteristic features of Proposition 1 to characteristics of the Mobius involu- 
tions, but can also be seen directly by considering formula (3.6) in the degenerate 
case when ri and r2 arc degree two polynomials. 

Expression (3.6) makes connection with Bezout's formula for the resultant of 
polynomials ri and r2 (see [27]), specifically the resultant is the determinant of 
the coefficient matrix of (3.6). Its appearance here is interesting, but perhaps 
not completely surprising given the role played by the resultant in Definition 1. 

As a final remark, we note that the pattern of terms and variables in (3.5) 
is very similar to that of the three-term sigma- function identity (see p. 390 of 
[28]), or Riemann relation, though an explanation for this is not so clear. 

The characterising polynomial identity (3.5) leads directly to the following 
constructive existence-and-uniqucness result for the idempotent biquadratics. 

Lemma 1. Given a pair of disjoint triplets taken from C U {oo} 

xo,xi,X2, 2/0,2/1,2/2, (3.10) 

there exists an idempotent biquadratic, which is unique up to a scalar multiple, 
for which these are orbits. It may be expressed via the formula (3.6) by choosing 

ri{x) = {x-xo){x-xi){x-X2), r2{x) = {x-yo){x-yi){x-y2). (3.11) 

Proof. For existence it is sufficient to check that the roots of the polynomials 
ri and r2 give two particular orbits of h in (3.6). Conversely, supposing (3.10) 
are the orbits of a given idempotent biquadratic polynomial h, choose w = xq 
and z — ya in the identity characterising such polynomials (3.5). By inspection 
it can be seen that the polynomials ri, r2 appearing in (3.6), (3.8), which is just 
a re-arrangement of (3.5), vanish on the orbits (3.10) due to the choice of w and 
z. In other words, any such biquadratic h coincides with (3.6), (3.11) up to a 
scalar multiple. 

Thus two orbits can be chosen freely, which then determine an idempotent bi- 
quadratic polynomial uniquely up to a scalar multiple. Clearly three orbits can- 
not be chosen freely. 

4. Rationality of the three-orbit constraint 

We connect the previously defined class of algebraic correspondences with the 
scalar quadrirational dynamical systems, this emerges directly by looking more 
closely at the three-orbit constraint. 

The following form of this constraint is a straightforward consequence of 
Lemma 1. 

Proposition 2. Denote three elements in the vector space of cubic polynomials 

by Ti} T2 and r3. The roots of these three polynomials are three distinct orbits of 
an idempotent biquadratic correspondence if and only if there are no pair-wise 
common roots and they are linearly dependent: 

371,72,73 G C\ {0} such that 7iri(a;) +'y2r2{x) -1-73^3(2;) = Mx. (4.1) 
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Fig. 5. Quadriratioiiality; variables on any pair of adjacent edges are determined rationally 
from the variables on the opposing edges. 

In terms of the roots of ri, r2, rs, i.e., directly in terms of the points of the three 
orbits, say 

Xq,Xx,x^, 2/0,2/1,2/2, zq,Zx,Zi, (4.2) 

the constraint leaves seven degrees of freedom amongst the nine variables. An 
important feature is its symmetry, by which wc mean invariancc under permu- 
tation of the three orbits, as well as permutation of the three points of any one 
of those orbits. The main reason for drawing attention to this symmetry is that, 
although transparent on the level of Proposition 2, it is obscured when writing 
the constraint in a more explicit form. 

One possible explicit form appears by speciiying all variables in (4.2) except 
two taken from the same orbit. For instance, z\ and are determined in terms 
of the remaining variables by a quadratic equation, hi^z^, z^ = z ^ {zi,Z2}, 
with h given by (3.6), (3.11). However, if two variables arc chosen from dis- 
tinct orbits, then they are determined rationally. To see this consider how the 
biquadratic (3.6), (3.11) depends on all data 

h{x,y;xo,xi,X2;yo,yi,y2) = 

{x - y)~'^[{x - xo){x - xi){x - X2){y - yo){y - yi){y - 2/2)- 

{y - xo){y - xi){y - X2){x - yo){x - yi){x - 2/2)]- (4.3) 

The dependence is quadratic in x and y, but in the remaining variables h has 
degree one. Writing the constraint on the three orbits (4.2) in the form 

h{xo,xi;zo,zi,Z2;yo,yi,y2) = 0, h{yo,yi; zq, zi, Z2;xo,xi,X2) = 0, (4.4) 

it is clear from the degree of the polynomials that 2/2 and X2 are determined 
rationally from the remaining variables. 

Consider now combining the rationality of the constraint with its permutation 
symmetries xi -H- X2 and yi ^ 2/2- It is clear that with xq, 2/0, zq, zi and Z2 fixed, 
the constraint on the remaining four variables xi, X2, 2/i and 2/2, visualised by 
assigning variables to edges of a quad as in Figure 5, determines any pair of 
adjacent variables rationally from those on the opposing edges. This system 
therefore falls within the quadrirational class introduced in [2]. An exhaustive 
list of five scalar quadrirational models was found, and denoted F/,...,Fy. 
The primary model, Fj, was written in the introduction (2.1), it is exactly the 
condition required for existence of an idempotent biquadratic polynomial with 
the three particular orbits 



u,u,/3, v,v,a, 0,1,00. 



(4.5) 
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Notice that the third orbit is constant. When two or all three points of the 
constant orbit coincide, the associated quadrirational model is equivalent to 
Fii or Fiv respectively. The remaining models, Fm and Fy, are obtained by 
limiting procedures that erode the connection to the idempotent biquadratics 
and are not considered here. 

The dynamics defined by the quadrirational models Fj, Fjj and Fjv are 
therefore associated with a family of idempotent correspondences that have a 
common fixed orbit, the distinction between the three models being the number 
of points in that orbit. This distinction is not important here, where possible we 
consider the generic form of these models (4.4) in which variables zq, z\ and zi 
of the third orbit are assumed fixed, but are left unspecified. 

The next step we take in this paper is the introduction of a second class of 
integrable dynamics, distinct from the quad-graph, which emerges from elemen- 
tary considerations of this family of algebraic correspondences. 

5. Oriented tetrahedral consistency 

We establish a closure/associativity property in the family of idempotent bi- 
quadratic correspondences that share a common orbit. This can be seen as a 
consistency feature of the rational three-orbit constraint (4.4) . It allows to in- 
terpret it as a discrete integrable dynamical system whose domain is the (edge-) 
oriented n-simplex. 

The relevant closure property may be stated as follows. 

Proposition 3. Denote twelve distinct elements taken from C U {oo} by 

Xij, ^,JG {1,2,3,4}, |{i,i}| = 2, (5.1) 

and take 

Xij , Xjkj Xki, XjiiXkjiXik, i, J, € {1, 2, 3, 4}, =3, (5-2) 

to he pairs of orbits that define ( of. Lemma 1 ) four idempotent biquadratic poly- 
nomials. If there exists an orbit common amongst any three of these, then it is 
common amongst all four. 

The geometric object we associate with this arrangement of variables is the 
tetrahedron with oriented edges. Indices are associated with vertices, variables 
with edges, and orbits with faces, see Figures 6(a) and 6(b). 

More algebraic intuition is offered by interpreting this property as associativ- 
ity. Consider the binary operation 

(a;o, t/o) • {xuVi) := (?/2, 2:2) (5.3) 

defined by system (4.4) when zq, Zi and Z2 are considered fixed. The commuta- 
tivity of this operation is clear from the permutation symmetry of the defining 
system, the associativity can verified directly by calculation. This is equivalent 
to verifying Proposition 3, which can be seen as follows. 

The proposition states that if the three-orbit constraint (4.4) is satisfied for 
three of the idempotents, then it is satisfied for the fourth. Taking the common 
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Fig. 6. (a) Tetrahedron, (b) one edge and two variables, one variable is assigned to each 
orientation of the edge, (c) a planar triangle tessellation. 

third orbit to be zq, zi, Z2, three of the constraints written in terms of the above 
binary operation are as follows 

(0:14,2:41) = (2:13, a;3i) • [xzi.Xiz], 

{x24, Xi2) = (a;2i, 2:12) • (a;i4, 2:41), (5.4) 
(2:23, 2:32) = (2:21, X12) • (2:13,2:31). 

Assuming these hold one sees by substitution that the fourth constraint, 

(2:24, 2:42) = (2:23, 2:32) • (2:34, 2:43), (5.5) 

is equivalent to the associativity condition 

(2:21,2:12) • [(2:13,2:31) • (2:34,2:43)] = [(2:21,2:12) • (2:13,2:31)] • (2:34,2:43). (5.6) 

That the equations represented by (5.5) are a consequence of those repre- 
sented by (5.4) means this is a consistency property similar in spirit to the known 
integrability feature of the quadrirational models, namely the consistency on a 
cube, or Yang-Baxter property. Here equations on three faces of the tetrahedron 
(5.4) imply the equation on the fourth face (5.5). 

The natural extension of this domain is the n-simplex. Coordinatization is 
immediate by replacing the range of indices {1, 2, 3, 4}, appearing the definition 
of variables (5.1) and orbits (5.2) of the 3-simplex, with the set {1, . . . ,n + 1}. 
Tessellation by triangles of any two-dimensional surface (see Figure 6(c)) can be 
considered as a sub-case of this domain, which is clear because the n-simplex 
can be recovered by adding edges until all pairs of vertices are joined. 

Models with such tetrahedral consistency can be directly integrated. Specif- 
ically the general solution of a given initial value problem can be written down 
in terms of the binary operation. Natural initial data for the oriented n-simplex 
domain is the set of pairs of variables associated with all edges connected to 
some vertex, for instance the n pairs (xij,Xji), j € {1, . . .n + l}\ {i} connected 
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20,21,^2 V^{ci,b) ip{a,b) ■ ip(c,d) 

0,1,00 (a(l - b)/{a - b), (1 - b)/{a - b)) ip{ac, bd) 

0,00,00 {a/{b - l),ab/{b - 1)) tp{a + c,bd) 

00,00,00 (a/b + b,a/b — b) ip(a + c,h + d) 

Table 1. Decomposing the binary operation into addition and multiplication. 



to vertex i. In this case the general solution in terms of the initial data is given 
by ixjk,Xkj) = {xj„x,j) ■ ixik,Xki) for any j,k G {1, . . . ,n + 1}\ {i}. 

The binary operation (5.3) defined by (4.4) is of itself very simple, it can be 
decomposed into addition and multiplication. Substitutions ip achieving this are 
listed in Table 1 for three canonical choices of the fixed third orbit. 

It is natural to ask what is the relationship between the dynamics observed 
here and the usual quad-graph dynamics. Answering this question is in fact the 
main technical focus of the remainder of this paper. We begin in the following 
section with a local consideration of this problem, establishing a combinatorial 
relationship between the cubic consistency and edge-oriented tetrahedral consis- 
tency properties. 

6. 5-simplex consistency 

Through a symmetrization procedure we introduce a particular domain for model 
(4.4) which is combinatorially equivalent to the 5-simplex. On this domain the 
usual cubic consistency (established in [2] for all quadrirational models) and the 
oriented-tetrahedron consistency (introduced in Section 5) are exhibited simul- 
taneously as sub-cases of a stronger consistency property. 

6.1. Completeness based on symmetry. It is useful to broaden the notion of do- 
main for (4.4) to a more combinatorial setting. Considered with the orbit zq, zi, 
Z2 fixed, the system (4.4) has the symmetry of a pair of disjoint triangles, there- 
fore in general the domain can be a set of variables arranged into triplet-pairs, 
where: 

Definition 2. A triplet-pair is a disjoint pair of sets, each set containing three 

variables. 

To remove degeneracy we assume connectedness, i.e., that any pair of variables 
in the domain are connected through some sequence of triplet-pairs. 

The domain appearing in Proposition 3 contains in total four triplet-pairs, 
whilst each variable participates in exactly two of them. A domain with this 
kind of regularity will be termed a pair- configuration. 

Definition 3. By a pair- configuration lue mean a set of variables arranged into 
triplet-pairs such that each variable participates in the same number of triplet- 
pairs. 

It means that throughout the domain, like within the equation itself, all vari- 
ables are on an equal footing. We remark that neglecting the pair associations 
between triplets recovers an abstract configuration in the usual sense. Also we 
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Fig. 7. The first diagram is the 5-simplex, fifteen variables are assigned to the fifteen edges, 
and ten triplet-pairs are assigned to the ten vertex-disjoint pairs of triangles, or polar triangles, 
a typical example of which is illustrated in the second diagram. 

deal exclusively with triplets here, but the more general definition would involve 
n-tuples of variables. 

By a symmetry of a pair-configuration we mean a permutation of the variables 
which also permutes the triplet-pairs. The group of symmetries of the simplest 
pair-configuration, a single triplet-pair (3.10), contains elements that permute 
each triplet, and also those that swap the two triplets. This group has order 
3! X 3! X 2 = 72. On the other hand, the oriented tetrahedral pair-configuration 
has just 2 X 4! = 48 symmetries, corresponding to the tetrahedral symmetry 
plus the overall reversal of orientation. By counting, it is therefore clear that not 
every symmetry of each participating triplet-pair extends to a symmetry of the 
full tetrahedral pair-configuration. For a pair-configuration which is a domain of 
our model we adopt the principle that this is a deficiency. To overcome it we seek 
to extend the configuration, such extension will be termed symmetry-complete. 

Definition 4. A pair- configuration is symmetry- complete when every symmetry 
of every triplet-pair is the restriction of a symmetry of the whole configuration. 

It is natural to require the extension be minimal, that is to ask for the smallest 
symmetry-complete pair-configuration that contains the original. The question 
of whether a finite completion of this kind exists is purely combinatorial. For 
the pair-configuration of Proposition 3 the answer is affirmative. 

Proposition 4. The smallest symmetry- complete extension of the edge-oriented 
tetrahedral pair- configuration (cf. Proposition 3) consists of fifteen variables 

Xij=Xju i,j e {1,2,3,4,5,6}, |{i,j}|= 2, (6.1) 

arranged into the ten triplet-pairs 

Xij,Xjh,Xki, XlmtXmnjXnU {^) J) ^) '^l ~ {!) 2, 3, 4, 5, 6}. (6-2) 

This extended configuration has the combinatorics of the 5-simplex, see Figure 
7. The indices are assigned to vertices of the figure, the variables to edges and 
the triplet-pairs to pairs of vertex-disjoint (or polar) triangles. 

Proof (of Proposition 4)- That the new configuration is symmetry-complete in 

the sense of Definition 4 can be seen from the figure, the symmetry group of 
the 5-simplex is 5*6 which acts naturally on its vertices. By inspection this group 
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acts transitively on the triangle-pairs, whilst the sub-group that stabilizes any 
particular triangle-pair is the complete symmetry group of that pair. 

To see that this configuration is minimal consider the intersection between two 
distinct triplet-pairs. In the case of the tctrahedral configuration the intersection 
is of a particular nature, triplet-pairs intersect on two variables - one from each 
triplet. For a hypothetical symmetry-complete extension of this configuration it 
is necessary that for any participating triplet-pair and any two variables chosen 
one from each of its triplets, that there exist another triplet-pair intersecting 
the first on those variables. For the 5-simplex pair-configuration described, a 
triplet-pair intersecting in this way exists and is in fact unique, and furthermore 
all intersections are of this kind. The extended configuration therefore cannot 
be any smaller, for if it were it would omit some of the necessary intersections. 

6.2. The stronger consistency property. The symmetrization procedure thus ex- 
tends the teterahedral domain to one which respects the symmetry of the model. 
It is a remarkable additional feature of (4.4) that the well-posedness of the ini- 
tial value problem is also preserved in this extension. In order to describe this 
situation precisely the following notion is useful. 

Definition 5. A duad-syntheme of the 5-simplex is a triplet of vertex-disjoint 
edges. 

This terminology follows Sylvester [29], a more precise usage would be edges 
associated with a duad-syntheme. The duads here are un-ordcrcd pairs from the 
set {1, . . . , 6}, and a duad-syntheme is a partition of this set into duads. In the 
same terminology a pair of vertex-disjoint triangles, to which we have previously 
assigned a triplet-pair, would be called a triad- syntheme. There are fifteen dis- 
tinct duad-synthemes in total, the associated variables are given explicitly as 
follows 

Xij, Xki, Xmn, {i,i,fc,/,m,n} = {1,2,3,4,5,6}. (6.3) 
The notion of syntheme allows concise statement of the consistency feature. 

Proposition 5. Consider model (4-4) on the 5-simplex domain (i.e., the domain 
described in Proposition 4)- Generic initial data on any pair of disjoint duad- 

synthem,es determines a unique solution. 

This can be verified directly by calculation. Only a single calculation is necessary 
because any pair of disjoint duad-synthemes can be transformed to any other 
pair by the action of the 5-simplex symmetry group. Note that this does not 
exhaust the initial data configurations that determine a unique solution. 

The consistency property may be understood more clearly through the follow- 
ing decomposition which again involves the duad-synthemes, but in a different 
manner. Distinguish any one of the duad-synthemes and partition the triplet- 
pairs of the 5-simplex pair-configuration into two sets according to whether or 
not they intersect it. There are four triplet-pairs disjoint from the distinguished 
syntheme, they lie in the oriented-tetrahedron pair-configuration of Proposition 
3. The remaining six triplet-pairs are in natural correspondence with the six 
faces of a cube with characteristics. To illustrate the association we have la- 
belled the edges of the 5-simplcx in Figure 8 according to the notation of the 
cube in Figure 3. The parameters, a, j3 and 7 lie on a duad-syntheme which 
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Fig. 8. The 5-simplex with variables assigned to edges. 



induces the following partition of triplet-pairs: 
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a, w 





(6.4) 



One; c;an verify the oriented tetrahedron pair-configuration on the left, whilst the 
six triplet-pairs on the right can be identified with faces of the cube domain of 
Figure 3. Consistency on the cube is usually described in terms of initial data 
Q:, /3,7 and u,v,w. This initial data allows to see equations on the on the left 
in (6.4) as a consequence of equations on the right. But there are fifteen ways 
to make such a partition, corresponding to the fifteen distinct duad-synthemes. 
There is therefore a harmony between the cube and tetrahedral consistencies, 
and it is natural to see the 5-simplex consistency as a stronger property of (4.4) 
which contains both as sub-cases. 



6.3. The multi- quadratic model Q4*. In this setting, which generalises the cube, 
the multi-quadratic model Q4* relates variables on the edges shown in Figure 9. 
This equation, returning to the notation of Proposition 4 (which reflects better 
the symmetry, corresponding to permutations of indices 1, . . . , 6 of the variables 
Xij), takes the following generic form 



xl){zo, Zi, Z2) + tp{zi,Z2,Zo) + ^p{z2,Zn,Zi) + 

ip{z2,Zi,Zo) +i>{zo,Z2,Zi) +1p{zi,Zo,Z2) = 0, (6.5) 
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Fig. 9. An edge of the 5-simplex and its polar tetraliedron. 

where 

1p{Zo,Zi,Z2) = iXi2 - Z2){Z0 - ) Cr(a;34 , Zi , X56 , ^2 ) cr(a;34, ^2 , 2^56 , -Zo) 

X (cr(a;45,xi2,.X36, Z2) - cr(x46, zi, .T35, zo)) (6.6) 
X (cr(3-;45, zi, 2:36, Zq) ^ cT:{xi(,, a;i2, X35, Z2)) , 

and cr denotes the cross-ratio 

cr(a,6,M) ■•= ^r'\^rl (6.7) 
[a — c)[o — d) 

It involves the variable X12, which is singled out, and variables on the edges 
of the tetrahedron formed by the remaining vertices 3,4,5,6. The tctrahcdral 
permutation symmetry of this model can be verified, though it is not obvious 
from (6.5), or from the expression (2.4), which corresponds to the particular 
case of (6.5) when X12 = c and {zq, 21,22} = {— c, 1/c, — 1/c}. The 5-simplex 
domain carries fifteen copies of this equation, corresponding to the fifteen edge- 
tetrahedron pairs, cf. Figure 9. 



6.4. A point-line configuration equivalent to the 5-simplex pair-configuration. 
We show a natural way in which the 5-simplex pair-configuration can be viewed 
as a point-line configuration in a three-dimensional projective space. 

Take the fifteen variables to be points and the set of all triplets from the 
triplet-pairs (6.2) and the duad-synthemes (6.3) to be lines (sets containing 
three points), so there are 35 lines in total. It can be verified that these are 
the 15 points and 35 lines of PG(3,2), i.e., the projective geometry in three 
dimensions defined over the field of characteristic two. The subset of 20 lines 
(6.2), whose union contains the entire set of 15 points, is the configuration of 
interest. Each point is on four lines and each line contains three points, it is a 
(I542O3) configuration. In PG(3, 2) each line is contained in three planes, and 
within this particular configuration there is a unique other line, skew to the first, 
which passes through all three of those planes. This relationship is reciprocal, and 
it recovers the pair associations of the pair-configuration. Whereas the system 
(4.4) lies on these pairs of skew lines, the equation (6.5) lies on planes. The 
complement of this configuration in PG{3, 2), i.e., the set of lines corresponding 
to the duad-synthemes (6.3), is known as the Cremona- Richmond configuration 
[30]. 
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7. The vertex-map group 

It is natural to ask what is the more general domain for models such as (4.4) 

which have the stronger 5-simplex consistency property. Formulation in terms of 
vertex maps allows group analysis to be applied to this question. To introduce 
the vertex maps we re- visit the oriented n-simplex and n-dimensional hypercube 
domains natural for models with the oriented-tetrahedral and cubic consistency 
properties respectively. 



7.1. Vertex-maps of the oriented n-simplex. The natural initial-value problem 
for system (4.4) on the oriented n-simplex domain (cf. Section 5) suggests the 
following elementary construction. 

Choose some vertex of the oriented n-simplex and denote the pairs of variables 
associated with its n connected edges by 



(7.1) 



the order in which edges are chosen is free, but we have chosen Xi to denote the 
departing variable and Ui to denote the arriving one (cf. Figure 6(b)). From such 
initial data the system (4.4) determines rationally all remaining edge variables 
of the n-simplcx, and in particular the data similar to (7.1) associated with any 
other vertex. In particular, on the new vertex connected to the initial vertex by 
initial edge j, this data is obtained by a mapping that we denote Tj, 



(7.2) 



where 

{xi,yi) = iyj,Xj) ■ {xi,yi), i e {1, . . . , n} \ {j}, 
= (yj^Xj). 

This mapping is expressed here in terms of the previously introduced binary 
operation (5.3) defined by (4.4). 

The mappings ti , . . . , t„ generate a group which we call the n-simplex vertex- 
map group. It is isomorphic to the full symmetry group of the n-simplex, A^. 
The generators satisfy the relations 





xi,yi 























(7.3) 



id, 
id. 



inTjTiTk)'^ = id. 



|{^,J}|=2, 
\{i,j,k}\ = 3, 



(7.4) 



which, it can be verified, constitutes a (non-standard) presentation for this group. 
The order of the group is (n -I- 1)!. If we consider the vertex data (7.1) only up 
to the ordering of connected edges (n! permutations), then the natural group 
action corresponds to permutations of the n -|- 1 vertices. 

The relations (7.4) encode the oricnted-triangle permutation symmetries and 
oriented-tetrahedral consistency property of (4.4). 
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7.2. Vertex-maps of the n- dimensional hypercube. The hypercube is the generic 
quad-graph domain for models with the cubic consistency property, containing 
for instance the planar quad-graph (with no self-intersecting characteristics) by 
a procedure of restriction. This domain also allows to define the vertex maps. 

Consider the following regular initial- value-problem. Denote by j/i, . . . , y„ the 
parameters associated with the n characteristics of the n-cube. Then select a 
vertex and denote by .xi, . . . , .t„ the variables on edges connected to that vertex, 
with ordering chosen so that the edge associated with variable Xj intersects the 
characteristic associated with parameter y^. Arrange the n variable-parameter 
pairs in the array (7.1), the overall ordering of pairs may be chosen freely. 

Prom such initial data the system (4.4) determines rationally all remaining 
edge variables, and in particular the similar data associated with any other vertex 
of the n-cube. On the vertex at the other end of initial edge j the similar data 
is obtained through a mapping which we denote by aj: 



(7.5) 



where determined in terms of polynomial h (4.3) by equations 

h{xj,yi]Zo,zi,Z2;xi,yj,Xi) = 0, i G {1, . . . ,n} \ {j}, 











Gj := 










Xni Vn 




Xnt Vn 



(7.6) 



The mappings cti , . . . , (t„ correspond to reflections normal to each axis of the 
n-cube, in particular they satisfy the relations 



o-i = id, 
(ai(7,)2 = id, |{i,i}|=2. 



(7.7) 



which encode the quad-symmetry and cubic consistency (or Yang-Baxter prop- 
erty) of (4.4). These mappings generate a group of order 2" and preserve the 
ordering of characteristics. Combined with pure permutation of the character- 
istics, i.e., the n! maps that permute the rows of (7.1), they generate the full 
symmetry group of the n-dimensional hypercube, BCn of order 2"n!. We call 
this the vertex-map group of the n-cube. 

It is possible here to make contact with the well-known braid-group formula- 
tion of the Yang-Baxter property. It corresponds to a similar construction to the 
one above, but a path-type initial value problem for the n-cube. Consider again 
specifying the parameters on characteristics but now initial data on a sequence 
of n edges connecting two (hyper-) opposite vertices of the n-cube. We write 
this data as an array of n edge-variables each paired with the parameter associ- 
ated with the characteristic that edge crosses, so the data takes the same form 
(7.1). The group of mappings that take one such initial data path to another are 
generated by those that take the simple form 



(7.8) 





xi,yi 
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where j G {1, . . . ,n — 1} and 

Xi = Xi, yi = yi, ie{l,...,n}\{j,j + l}, 

h{xj+i,yj+i;zo,zi,Z2;Xj,yj,Xj) = 0, y, = yj+i, 
h{xj , yj ■,zo,zi,Z2;Xj+i, y^+i ,Xj+i) =0, yj+i = yj . 

They correspond to deforming the path by flipping edges across a single quad, 
in other words interchanging adjacent characteristics. The characteristics cor- 
respond to braids and the Yang-Baxter property is encoded in the first of the 
braid relations 

SiSjSi = SjSiSj^ \i j\ = 1, 

SiSj = SjSi, \i - j\ > 1. 

Here sf — id, so for these systems the braid relations lead only to the symmet- 
ric group. The intended extension beyond this group (which means extension 
beyond the hypercube) is initially less transparent in terms of the braid-type 
mappings. On the other hand, the hypercube geometry makes it straightforward 
to construct a transformation p that takes the vertex-type initial data to the 
path-type initial data. Denote by (j) the column vector (7.1) and define p by the 
expression 

p-- (j)^ diagonal([(/),CTi(0), [cricr2](0), • • • , [tri • ■ ■ an-i]{(l>)]). (7.11) 

The inverse is obtained similarly, 

p-\(b) - first-row([0, .si((/)), [siS2](0), . . . , [si • • • .s„_i]((/))])^. (7.12) 

Note that arrays in (7.11) and (7.12) are to be understood as n x n square 
matrices whose entries are pairs of variables such as {xi,yi) etc. Prom this 
explicit form for p it is possible to verify, for instance, that the transformed 
braid-type maps, 

p~^Sip, p~^Sn-lp, (7.13) 

are standard generators of the pure permutation part of the hypercube vertex- 
map group. 

Note that from here it is possible to make contact also with the transfer 

dynamics of Veselovs monodromy maps [31]. In terms of the lattice geometry it 
corresponds to tessellation of Z" by hypercubes, and the monodromy maps are 
identified with dynamics of the periodic, or finite-gap, reduction (details of this 
were established in the context of more general intertwining Yang-Baxter maps 
by Kouloukas and Papageorgiou [32]). 

7.3. Characterisation of the complete vertex-map group. The previous two sub- 
sections establish transformation groups that act naturally on vertex-type initial 
data of the n-dimensional edge-oriented simplex and hypercube respectively. 
Combining both sets of generators in the case n = 3 yicilds a group of trans- 
formations (a-i,0-2,(T3,Ti,r2,r3) which, it can be verified, is the transformation 
group for initial data configurations of the 5-simplex pair-configuration of Propo- 
sition 5. In this case the variables xi, X2, .xs and yi, y2, ys of (7.1) are associated 
with a pair of disjoint duad-synthemes (cf. Definition 5) of the 5-simplex (despite 



(7.9) 



(7.10) 
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Fig. 10. Coxeter-Dynkin diagram for the index-2 normal subgroup of the vertex-map group. 



this we retain the phrase vertex maps, which recalls the original context of the 
construction). Thus combining the vertex maps provides a constructive method 
to recover the 5-simplex domain that was obtained in Section 6. On this basis, 
we now investigate the group generated by the combined vertex maps for generic 
n. 

Analysis of this group is simplified by the observation that vertex maps from 
the two classes can be related, 

Tj=ajuaj, (7.14) 
by introducing the map w that interchanges columns 





' xi,yi 




yi,xi 


CO := 










^n; Vn 







We assume the integer n is fixed and consider the group generated by 

ai,...,am,uj, (7-16) 

for m < n. This has the advantage that if m' < m, then (cri, . . . , (7^', is a 
subgroup of ((Ti, . . . ,crm,,w). The following relations are key: 

af =ijp- = [oiiS)'^ = id, 
(Oicj.f = (a,c.a,w)3 = id, |{z, j}| = 2, (7.17) 
{{piajLofak^Sf = id, j, A;}| = 3, 

they encode the stronger 5-simplcx consistency property. 

The main result of this section is the characterisation of this group for generic 
m < n. 

Proposition 6. The finitely presented group with generators (7.16) and rela- 
tions (7.17) contains an index-2 normal subgroup isomorphic to the Coxeter 
group of Figure 10. 

Proof. The index-2 normal subgroup can be written as 



(w, cTiwcTi, cricr2, CTiCTs, . . . ,a\am)- 



(7.18) 
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Group Order Isomorphic to 

((Ti, 0-2,0; > 2 X 36 Ai K {A2 X A2) 

(cti , (72 , 0-3 , tj) 2 x 720 Ai K As 

(cri,(T2,(T3,(T4,u;) 2 x 51840 Ai tK Eq 

{cri,o-2, 0-3, <T4, 0-5,0;) 2 X 00 Ai x Ee 

Table 2. Sequence of subgroups of the complete vertex-map group. 



The claim is that this group is isomorphic to the group represented in Figure 
10. We proceed by exhibiting an isomorphism: 





= w, 


LO = ti, 






= aicjai, 


aiujai ■ 


= t2, 


h 




(J\a2 = 


tat 4, 


u 


= (wcri(T2)^, 


(Tiers = 




h 




criC74 = 





(7.19) 



tm+2 = {0Jam-l(Tmf, OxCfm = (i3t4)*'*'''"*'"+' • 

Here the notation denotes conjugation of / by 5, P — g~^fg. The group 
elements ti,. . . ,tm+2 arc associated -with the labelled nodes of Figure 10. An 
advantage of the chosen isomorphism is that to verify it, calculation is required 
only up to the case m = 4, because beyond this the new relations which appear 
are similar to existing ones that have already been checked. This can be seen 
by inspection of the isomorphism (7.19), the relations (7.17) and the Figure 
10. At the value m = 4 the groups in question are finite and the isomorphism 
may be checked using computation based on a permutation representation. The 
computer algebra system Magma [33] was used to construct an isomorphism at 
m = 4. 

Note that mappings of the form (wajaj+i)^, which appear in (7.19), correspond 
to pure permutations in the vertex-map group, specifically transposition of rows 
i and i + 1 of the array (7.1). In particular, we note that for i G {1, . . . , m — 1}, 
= (ujaiam-i)'^ = s^, where si, . . . , s„i_i are the braid-type mappings. 
We also remark on one subtlety, namely that m < n is required for the 
relations (7.17) to be complete. In fact, due to the groups being finite, it is easy 
to check that the following holds. 

Lemma 2. If m = n G {3,4} then (cti, . . . ,cr„, w) = {co, aicoai,aia2, ■ ■ ■ , cri0-„). 

This would contradict Proposition 6 if the relations were complete in the case 
TO = n. It may continue to hold for all to = n > 2 (but is not true if to = n = 2). 

The group relations (7.17) and Proposition 6 can be seen as a fairly complete 
characterisation of the algebraic structure associated with the family of idem- 
potent biquadratics that share a common orbit. To understand what this means 
for the associated discrete dynamics, the group structure requires some further 
interpretation. Within Figure 10 one can see the symmetry group of a single 
triplet-pair corresponding to to = 2, and of the 5-simplex pair-configuration 
when TO = 3. Continuing, the combined vertex-map group decomposes as in Ta- 
ble 2. In particular, that the group becomes infinite beyond to = 4 represents 
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a significant and un-anticipated depth beyond the kind of lattice geometry of 
the oriented n-simplex or the n-dimensional hypercube. A natural step towards 
understanding the dynamics of the ambient space as a whole, is to relate the 
lattice geometry to the group. 

8. The generalised lattice geometry 

The abstract structure of the vertex-map group encodes the full ambient do- 
main on which it is natural to consider model (4.4). Here we explain how the 
geometry of the domain may be decoded from the group, make contact with the 
configuration of the 27 lines of a cubic surface in the case n = 4, and discuss the 
dynamics associated with subgroups. 



8.1. Decoding the geometry from the group. Recall from Section 6 that in general 
the domain of (4.4) is a set of variables arranged into triplet-pairs (cf. Definition 
2). Consider now the following associations, 

variables <H- F := {af : g G (ai, . . . ,(7„, w)}, 
triplet-pairs -H- T := {(cti, (72, w)^ : 5 € (cti, . . . , a-„, w)}. 

Again here denotes conjugation of / by g, = g~^fg- Thus a set of variables 
V and triplet-pairs T are identified with particular conjugacy classes of the vertex 
map group. It may be decided whether a variable is contained within a triplet- 
pair in the natural way: the triplet-pair corresponding to tt e T contains the 
variable corresponding to a & V when cr € tt. In particular, the internal structure 
of the triplet-pair corresponding to n G T may be recovered by taking the 
intersection V (In and considering how the action of n permutes these elements 
of V. 

Note that all of cti, . . . , (j„ are conjugate (this can be seen from the definition 
of aj (7.5), (7.6) by considering the pure-permutation sub-group), therefore all 
are identified with variables. The vertex-type initial data can also be recovered 
from the group, it corresponds to variables associated with the set 

{ai,...,(Tn,(Ti,...,a'^}, (8.2) 

or any set conjugate to this one. Furthermore, arranging this set in the form 



.0-„,cr„ 

provides a method to see the consistency of identifications (8.1). One checks that 
the natural action of the group on (8.3) (i.e., by conjugation of each component), 
is consistent with its default action on (7.1). It suffices to check the action of the 
generators. 

The combinatorial considerations raised in Section 6 are also taken care of 
automatically by this construction. 

Proposition 7. The variables and triplet-pairs defined by associations (8.1) 
form a symmetry- complete pair- configuration (cf. Definitions 3 and 4)- 



(8.3) 
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Proof. The domain is clearly a pair-configuration because all members of V and 
T are on an equal footing amongst themselves (are conjugate) . Now observe that 
the symmetry group of the triplet-pair associated with w G T can be identified 
with the inner automorphism group of tt. Clearly the action of tt extends natu- 
rally to the whole group, and thus the pair-configuration is symmetry-complete. 



8.2. Case n = 4. The case n = 4 is significant, being the last finite case in 
the sequence of complete vertex-map groups, cf. Table 2. This group, which is 
isomorphic to Eq when n = m = 4 (Table 2 shows n > m = 4, cf. Lemma 2), 
was studied by Jorden [34] as the group of incidence-preserving permutations 
of the 27 straight lines on a cubic surface, which he connected with the trisec- 
tion of periods of quadruply-periodic functions. Subsequently Burkhardt [35], 
through the theta-function approach introduced by Jordan, obtained a linear 
representation of the group, and in particular an invariant of the representation 
that, when set equal to zero, defined a quartic hypersurfacc (quartic primal). 
Amongst the substantial literature on this subject, we cite here mainly the work 
of Baker [36] who studied the geometry of Burkhardt's quartic. This study was 
purely in terms of position in four-dimensional space of the 45 singular points, 
called nodes, which he found were in natural correspondence with 45 particular 
elements of the group, called projections. The set of projections corresponds in 
turn to the set V in (8.1), i.e., the variables of the n = 4 pair-configuration. The 
projections were proven to form a complete conjugacy class by Todd [37], who 
also makes a concise summary of Bakers extensive geometric description of the 
45 nodes. 

We recount a part of this geometry in order to identify in Bakers scheme 
what we have called triplet-pairs here. The 45 nodes lie by threes on lines, called 
K-lines, with 16 lines through each node and 240 lines in total. The nodes and k- 
lines therefore form a (45i62403) configuration. The set of all K-lines intersecting 
a particular one contain every node of the configuration except three. Those 
three remaining nodes themselves lie on a second K-line, which is called polar to 
the first. The relationship is reciprocal, and the 120 pairs of polar K-lines can be 
identified with the triplet-pairs T in (8.1). 

Within the n = 4 pair-configuration there are 36 distinct copies of the n = 3 
(5-simplex) pair-configuration. This sub-configuration is left unnamed by Baker, 
but is called a 15-set by Todd. On the other hand, one of Bakers additional 
contributions was to provide a coordinatization, and this implicitly acknowledges 
the 15-set because it directly extends the natural coordinatization of it given 
earlier in Proposition 4. It allows to write down the n = 4 pair-configuration as 
follows. 

Proposition 8 (notation due to Baker [36]). The pair-configuration defined 
by associations (8.1) in the case n = 4 consists of 45 variables arranged into 120 
triplet-pairs. If we denote 15 of the variables by 



i,ie {1,2,3,4,5, 



6},|{i,i}|=2, 



(8.4) 



and 30 of them by 



^{ij,kl,mn) ^{ji,kl,mn) ^{kl,mn,ij) ■> J5 ^5 ^} {l?^!^;'^?^? 



6}, (8.5) 
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then the 120 triplet-pairs may he written as follows: there are 10 of the form 

X^j ^ Xjf^^ Xfj;i^ Xijji^XjYinf'^nli (8*6) 

90 of the form 

and 20 of the form 

'^{ij,kl^7nn) i '^{il.kn^mj) t ^(in^kj^ml) 7 [ij ^mn.kl) t •^(il.mj.kn) 7 '^{in.nil.kj) 1 (8-8) 

where in (8.6), (8.7) and (8.8), {i,j,k,l,m,n} — {1,2,3,4,5,6}. The following 
set of variables may be given as initial data for this pair- configuration 



(8.9) 



Xl2 , 2^54 
X34 , XiG 

-^(12,34,56) ; 2^(16,54,32). 

The configuration also carries in total 540 = 36 x 15 copies of model (6.5): 
the pairwise intersection of the 36 copies of the 5-simplex sub-configuration is 
of two kinds, cither a triplet-pair, or a duad-synthcme, thus no two share the 
same edge-tetrahedron pair of Figure 9. We also remark there are in total 135 
distinct copies of the 4-cube, which is the quad-graph domain generalised by this 
configuration. 



8.3. Subgroup dynamics and integrability. By construction, the complete vertex- 
map group contains An and BCn subgroups corresponding to the n-simplex and 
n-cube domains. In both of these domains it is envisaged that n is arbitrarily 
large, and that the number of degrees of freedom is therefore un-bounded. In- 
tegrability on the quad-graph domain is of KdV-type, whereas the integrability 
on the triangle-graph domain (discussed in Section 5) is of a simpler kind. 

Due to the nature of the complete vertex map group, the unbounded domain 
required for non-trivial dynamics is also present for finite n, the simplest case 
being n = 5. The group in this case is the afRne reflection group Eq, such afflne 
groups arc familiar from dynamics in the Painlevc setting, though the birational 
representation here looks very different. Nevertheless it would be natural to 
adopt the similar point of view, that dynamics are associated with the normal 
abelian subgroup, whilst the action of the finite quotient group (Eq) corresponds 
to symmetries of the dynamical system. This would give a natural notion of 
integrability, but it remains open whether the systems here can actually be 
integrated in terms of (discrete) Painleve transcendents or their generalisations. 

Restriction to subgroup dynamics by taking a reduced set of generators leads 
to an additional variety of possibilities. An interesting case would be the group 

{ti,t3,...,tio), (8.10) 
written in terms of elements given in (7.19), which constitutes a birational repre- 
sentation of Eg. Note that when restricting to a proper subgroup the associations 
(8.1) will in general require modification. 

Due to the rich structure of the general vertex map group, it is natural to 
expect that subgroup dynamics more generally can also be considered integrable 
in some sense. This is something which could potentially be tested by the notion 
of algebraic entropy [38,39] for instance. 
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9. Discussion 

Non-metric domains. It is relevant here to address the fact that the ambient 
domain, a symmetry-complete pair-configuration (Definitions 3 and 4) a-priori 
carries no natural notion of distance, wc loosely refer to such a domain as being 
non-metric. To use the language of dynamical systems can begin to be awkward 
in such situations. On the other hand, this domain is encountered here as the 
natural destination along the path of the consistency framework for integrable 
systems. What is most striking is the inherent potential of such domains to unify 
integrable systems from quite different settings. In particular, the one model here 
encompasses dynamics associated with a quad-graph, a triangle-graph, afRne 
reflection groups and more (cf. Section 8.3). 



Symmetry- com,pleteness and consistency. Here we have identified the consis- 
tency configuration underlying the integrability of two quad-graph models as 
being in fact a sub-configuration of a larger whole. Aside from the algebraic 
motivation for identifying this generalised picture, which is based on the idem- 
potent biquadratic and therefore specific to the particular models at hand, it 
is shown that it can also be discovered through the perfectly transferable and 
wholly combinatorial concept of symmetrization (cf. Definition 4 and Proposi- 
tion 4). Moreover, the technique based on vertex maps (more generally, mappings 
between initial data sets of the configuration) provides a method to obtain the 
symmetrized configuration. These methods consider the consistency of itself and 
therefore have application beyond the system studied here. 



Geometry. Adler, Bobenko and Suris. when introducing the scalar quadrira- 
tional models [2] , gave also a geometric incidence picture in terms of pencils of 
conies. The refinement of this geometric picture is potentially very interesting. In 
particular offering the possibility to make a connection with the configuration of 
27 lines on the cubic surface beyond the purely combinatorial relationship given. 



Generalisations. Underlying the system studied here are the idempotent bi- 
quadratics. More specifically, the family that share a common fixed orbit, that 
we have denoted Zq,zi,Z2 € C U {oo} throughout. The invariance under per- 
mutation of orbits raises the question of whether there is a yet more general 
consistency configuration in which this fixed orbit participates on the same level 
as all the others. 

Beyond the idempotent biquadratic is the idempotent bicubic. This is defined 
by the similar Bezout-type expression (3.6), but with n and r2 now quartic 
polynomials. As a correspondence, it directly generalises the biquadratic case, 
and the three-orbit constraint is again rational. This is therefore a potentially 
significant object and will be the subject of future work. 
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A. Idempotent biquadratic polynomials in the Weierstrass 
theory of elliptic functions 

Consider the symmetric biquadratic polynomial, 

h{x, y) = co + ci{x + y) + C2xy + 03(3;^ + y^) + C4xy{x + y)+ c^x^y'^, (A.l) 

co,...,C5 G C. As mentioned in the main text, its orbits (3.2), (3.3) may be 
taken as a definition of elliptic functions in the discrete setting. We develop here 
an aspect of this discrete theory connected to the Weierstrass addition formula, 
and relate it to the idempotent class of biquadratics. 

A basic role in this aspect of the Weierstrass theory is played by the discrim- 
inant polynomial r associated with h, which is defined as 

/^[h{x,y),y\ = r{x), A := discriminant. (A-2) 

This terminology is unambiguous due to the symmetry of h. Henceforth we 
assume that r has at least one simple root. The fundamental property on the 
level of biquadratics is a kind of closure under composition within the family 

with shared discriminant polynomial. More precisely, if biquadratics hi and ^2 
have, up to a scalar multiple, a common discriminant polynomial r, then the 
composition of hi and /12 takes the form 

ves\h\{x, y), h^iy, z), y] = h3{x, z)h4{x, z), res := resultant, (A-3) 

where /13 and /14 are again symmetric biquadratic polynomials whose discrimi- 
nant polynomial is a (now possibly zero) scalar multiple of r. The appearance 
of two factors on the right-hand-side of (A. 3) means this composition gives a 
structure which is sometimes called a two-valued group. 

A natural parameterisation of this group is obtained through the evaluation 
map, 

/i H> a : ft,(eo, a) = 0. (A.4) 

Choosing the point of evaluation, eo € C U {00}, to be the assumed simple 
root of r, makes the evaluation map single-valued. Furthermore, h can be re- 
constructed (up to constant scalar multiple) from its discriminant polynomial 
and the associated parameter a. Specifically 

h{x,y) ^ H{x,y,a) (A. 5) 

where if is a triquadratic polynomial which may be expressed as a discriminant 

t- eo J 

where f is a dummy variable, and 

n{t)^p^, r2{t) = {t-x){t-y)(t~a). (A.7) 

Observe that H is invariant under permutations of the three variables x, y and 
a. This permutation symmetry has the remarkable consequence that the param- 
eters ai, . . . , a4 associated with the biquadratics appearing in the composition 

formula (A. 3), 

hii-^ ai : hi{eQ,ai) =0 ^ hi{x,y) H{x,y,ai), ie {1,2,3,4}, (A.8) 
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are themselves related by the same polynomial H, 

H{ai,a2,z) = <^ z€{ai,ai}. (A.9) 

In particular, the binary operation (x, y) z defined by equation H{x, y,z) =0 
endows C U {oo} with the structure of an abclian two-valued group. Most of 
the important features of H (A. 6) can be seen from its canonical form, the case 
eo = oo, r{t) = At^ - g2t - gr- 

H{x, y, z) = {xy + yz + zx + \g2f - 4(x + y + z){xyz - igg). (A.IO) 

Those features being that it appears as the discriminant of a polynomial which is 
degree-one in each of three variables, that it is symmetric in the three variables, 
and, through the Weierstrass addition formula 

if(p(a),p(/3),p(7)) = ^ p(7) = p(a±^), (A.ll) 

also that it defines a two- valued abclian group structure. The more generic form 
(A. 6) reveals the connection with the Bezout-type expression for the idempotent 
biquadratic correspondences (3.6). 
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